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ABSTRACT 

We study the consequences of a homogeneous dark energy fluid having a non-vanishing velocity 
with respect to the matter and radiation large-scale rest frames. We consider homogeneous 
anisotropic cosmological models with four fluids (baryons, radiation, dark matter and dark 
energy) whose velocities can differ from each other. Performing a perturbative calculation up 
to second order in the velocities, we obtain the contribution of the anisotropies generated by 
the fluids motion to the CMB quadrupole and compare with observations. We also consider 
the exact problem for arbitrary velocities and solve the corresponding equations numerically for 
different dark energy models. We find that models whose equation of state is initially stiffer than 
radiation, as for instance some tracking models, are unstable against velocity perturbations, 
thus spoiling the late-time predictions for the energy densities. In the case of scaling models, 
the contributions to the quadrupole can be non- negligible for a wide range of initial conditions. 
We also consider fluids moving at the speed of light (null fluids) with positive energy and show 
that, without assuming any particular equation of state, they generically act as a cosmological 
constant at late times. We find the parameter region for which the models considered could be 
compatible with the measured (low) quadrupole. 
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1 Introduction 



The present observational evidence from supernovae type la [1, 2], CMB anisotropies [3] and 
large scale structure, mainly through the baryon acoustic oscillations, suggests that today the 
universe could be dominated by a negative pressure fluid [4, 5]. Although these data are 
compatible with the presence of a cosmological constant, the fact that such observations only 
explore relatively recent epochs implies that other possible models in which the equation of 
state of dark energy could have changed in time cannot be excluded a priori. Thus, if dark 
energy can be parametrized as a perfect fluid with equation of state pde — W]:,e{z)pde, the 
above mentioned observations only constrain the present value as — — 0.97lo!o9 [3]- The 
redshift dependence of the equation of state can be parametrized in different ways and still a 
wide range of variability is compatible with observations [6]. Moreover, apart from including 
new components in the energy-momentum tensor within Einstein gravity, or new scalar fields 
as in the quintessence models [7], it has been also suggested that the accelerated expansion of 
the universe could be due to modifications of the gravitational action at large distances [8] . 

The possibility of finding observational signals of dark energy which could discriminate 
between the various models is thus becoming of crucial importance. However, up to date, only 
a few proposals have been considered in the literature, which can be broadly classified in two 
classes [9]: on one hand probes of the redshift-distance relation, as for instance the already 
mentioned high-redshift supernovae la or the use of baryon acoustic oscillations as standard 
rulers. On the other hand we have the probes of the growth of structure in the universe, such 
as the weak gravitational lensing or the integrated Sachs- Wolfe effect. The goals of the future 
observational surveys will be the determination of the equation of state of dark energy with a 
few percent accuracy and the possibility of discriminating from a pure cosmological constant 
[10]. 

In this paper we consider a different aspect of dark energy with potential observational 
consequences which is related to the possible motion of dark energy with respect to matter and 
radiation. If dark energy can be described as a homogeneous perfect fiuid then, apart from 
the density parameter Q,de and equation of state wde, a complete knowledge of its energy- 
momentum tensor requires the determination of its relative velocity with respect to the rest of 
components of the universe. Indeed, dark energy is usually considered as a highly homogeneous 
fiuid in a similar way to radiation. This is due to the fact that in most models, its sound speed is 
close to the speed of light and this prevents the growth of dark energy perturbations below the 
Hubble scale. Moreover, dark energy is also required to be extremely weakly interacting with 
baryons and radiation, in order to avoid confiicts with the predictions of Standard Cosmology. 
Indeed, in most of the models, its effects are purely gravitational and dark energy is considered 
as a totally decoupled fiuid. In such a case it makes sense to ask whether the dark energy rest 
frame converges towards the radiation or matter large-scale rest frames. Since there is no a 
priori reason to expect that dark energy were necessarily coupled to radiation in the very early 
universe, its initial velocity with respect to radiation should be considered as a free cosmological 
parameter on equal footing to Qde or Wde, and very much in the same way as those parameters, 
the relative velocity of dark energy should be determined by observations. In addition, the fact 
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that a pure cosmological constant is invariant under change of frame imphes that the potential 
effects associated to a non-vanishing relative velocity will be exclusively present in models with 
varying equation of state. 

The metric anisotropies generated by the fluids motion can affect both the temperature and 
polarization of the CMB. In a previous work [11], we have started the study of such effects 
on the CMB dipole. We have shown that the motion of dark energy is not incompatible 
with the current measurements of the dipole, but instead, it modifies its usual interpretation. 
Thus, when dark energy is moving, the dipole is generated by the motion of emitter and 
observer with respect to the cosmic center of mass and not with respect to the background 
radiation. This fact can also have important consequences for the generation of matter bulk 
flows on very large scales [12]. In this work we consider the effects of motion on the CMB 
quadrupole. In the last years, the interest in anisotropic models (see [13, 14] and references 
therein) has grown mainly motivated by some unexpected features in the low multipoles of 
the CMB temperature anisotropies, in particular the low value of the quadrupole and the 
quadrupole-octupole alignment which could suggest the existence of a preferred direction in 
the universe [15, 16]. In this work we explore this possibility by considering cosmologies with 
several moving fluids (see [17] for previous works). 

The paper is organized as follows: in Section 2 we obtain the metric solution for a model with 
moving fluids up to second order in the velocities. In Section 3, we obtain the corresponding 
CMB temperature anisotropies and the contribution to the quadrupole. In Section 4 we consider 
the exact problem for large velocities and obtain the general expression for the quadrupole. 
Section 5 is devoted to the application of the previous results to different dark energy models. 
Finally Section 6 contains the main conclusions of the paper. 



2 Slow-moving fluids: second order equations 

Let us consider a universe filled with four homogeneous perfect fiuids: baryons, radiation, dark 
matter and dark energy. For matter and radiation we shall consider their usual equations 
of state, i.e., pr = and pb = Pdm = 0, whereas for dark energy we shall use: pde ~ 
wde{z)pde-i which in general depends on redshift. In an arbitrary frame, the energy-momentum 
tensor for each component takes the perfect fluid form: 

T!^''^{pa-rPo)<K-p^g^\ (1) 

where a — B,R, DM, DE. We shall consider homogeneity in the fluids so that all the quantities 
appearing in (1) will just depend on conformal time: 

Pa = Pair]), 

< = la{v){l,Uv)) (2) 

with 

la = I ^ . (3) 

Vfl'oo + Qij^lvi 
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In order to simplify the problem and to obtain some analytical solutions we shall use perturba- 
tion theory, assuming that the fluids velocities are small, i.e. <^ 1. To that end, we expand 
the different quantities of the four fluids up to second order as follows: 

Pa = PL°)+P«+Pi^) + --- 

V. = ^^^+^^^ + --- (4) 

where we have imposed the fluids to be at rest to zeroth order, i.e., vj^^ = 0. That way, the most 
general form for the metric is given by the perturbed Friedmann-Robcrtson- Walker metric: 



ds — a 



1 + 2 (0(1) + dri^ + 2 [^f ) + S^] dx'dri- 
(l - 2 (V'(i) + V^^^)) + hij\ dx'dx^] (5) 

where we will follow the notation in [18]. In this expression, (p^^^ and t/''-^-* are scalar pertur- 
bations of flrst order and will be determined from p^^^ in the flrst order equations of motion. 
However, the second order scalar perturbations 0^^^ and ip^"^^ can depend, not only on p'^^ and 
{Pa^y terms, but also on {v^^^'Y which are also scalars. Analogously, the flrst order vector 
perturbations S^^^ can only be related to the first order velocities v^'^ in the equations of mo- 
tion, whereas to second order, 5^*^^^ will be determined by combinations of v^^ and F^^^v^'^ 
where F^^^ is a scalar function of the first order scalar perturbations. Finally, hij is a trace- 
less tensor perturbation which should be of second order and depend on the combinations 
'^ai'^aj ~ ^i'^a^y^ij- Note that, slncc we are considering only time-dependent perturbations on 
the energy-momentum tensor, all the perturbations on the metric will be functions only of time 
and, therefore, the perturbed metric does not contain any terms involving spatial derivatives. 

So far, we have not done any specific gauge choice, so we still have four gauge degrees of 
freedom. Hence, we can simplify the problem by chosing our coordinates appropriately. In 
particular, we can fix the spatial coordinates so that the vector part of the metric vanishes 
S — 0. The physical interpretation of this condition is apparent when solving the (° J Einstein 
equation of the exact problem. Thus, we obtain the condition: 

^ ^ j:all{Pa+Pa)gijVi 
J2a7l{pa+Pa) 

where the metric is ds'^ = g^^dx^^dx^ . Hence, S can be interpreted as the relativistic cosmic 
center of mass velocity (see [11]). Notice that in general, an observer at rest with respect to 
cosmic center of mass could be moving with respect to radiation or matter. On the other hand, 
the temporal coordinate can be chosen in such a way that Z]a(Pa ~ P^a^) — 0, which means that 
the total density perturbations are identically zero. With this gauge choice, the q) and j) 
components of Einstein equations G^^^ = SttGT^j^ up to second order adopt the form: 



Zeroth order 



- ^^a^Y.p'^ (7) 

a 
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First order 



a 



(9) 
(10) 



Second order 



2n (V'(') + (V^'('))')' - (0(2) - 2 (0(1))') + ^'(1) (V^'(i) + 47Y0(i)) 



(11) 



(27i' + 7^2) 0(2) + ^^/(2) + ^//(2) ^ 27iV''(2) - 2 (7^2 + 27i') (0(i))' 



+ i^'(i) (^(1) _ 20(1))' + 2<(i) (^(1) - 0(1)) + 2n 
+ ^ {-'h'))' = SttG E (p(°) + p(°)) .(D'.W 



(i))2 + (0(i))2y + 20(^(1) 



(12) 
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with ' = and 7i — a'/a is the Hubble parameter. 

Nevertheless, the system is incomplete because there are more unknown variables than 
equations. In general, the problem with n fluids has ten independent Einstein equations, but 
the unknown quantities are the densities (assuming a given equation of state) and the three 
independent components of the four- velocity of each fluid (because of the constraint u'^ = 1). 
Therefore, there are 10 + An unknown functions, although, since there are four gauge degrees 
of freedom, we can fix four quantities and reduce the number of undetermined functions to 
6 + 4n. With this count, one needs 4(n — 1) additional equations to complete the system. The 
simplest way to close the problem is by requiring the conservation of each energy-momentum 
tensor, assuming they are decoupled from each other: T^^\ = 0. Obviously, one can modify 
these relations by changing the right hand side in order to consider interactions between the 
fiuids. This guarantees the completeness of the system since it provides the 4(n — 1) required 
equations (there are 4n extra equations, but the conservation of the total energy-momentum 
tensor makes one of those equations superfluous). For our case, these additional equations read 
for the energy and momentum conservation (notice that momentum conservation is trivial at 
zeroth order): 



Zeroth order 



First order 



pr+37i(p(°)+p(°)) = 
p'i^^ + 37i(pW+pW) = 3(p(°)+p(°))^'(i) 



(13) 

(14) 
(15) 
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Second order 



(2) , J2)^ 



+ 3 



(16) 



where we have used the previous orders equations at each order. 



3 Contributions to the CMB quadrupole 



The relevant part of the metric for the quadrupole is given in our case just by hij, since it 
is the only part contributing to the anisotropy. The homogeneous scalar perturbations only 
affect the value of the monopole in an negligible way. Concerning vector perturbations, we 
are working in the cosmic center of mass frame and accordingly those contributions vanish in 
our calculations. As commented in the introduction, the S contributions have been studied in 
[11] and modify the usual dipole contribution. Therefore, in order to calculate the quadrupole 
produced by this metric, we can consider just the tensor perturbation. Then, from now on, we 
shall use the metric: 



ds^ = c^iydrf — {Sij + hij)dx^dx^). 



In order to calculate all the contributions to the temperature anisotropics generated by the 
metric perturbations, we should solve the corresponding radiative transfer equations (see [19, 
20]). This is the system of Einstein-Boltzmann equations for the set of fluids. However, since 
we are only interested in the quadrupole anisotropy (which is not affected by microphysics at 
the time of recombination), the only relevant contribution for such a large angle contribution 
would be given by the Sachs- Wolfe effect which takes into account the variation in the energy 
of photons propagating from the last scattering surface [20]: 



6T 



(19) 



dec 



Here, the indices and dec denote the present and decoupling times respectively and E is the 
energy of the photon. For an observer with velocity — 7(1, v) this energy is given by: 



with 



dx" 



(20) 



(21) 
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where E parametrizes the photon energy and A is an affine parameter. By the invariance of 
the action of the geodesies of a massless particle under conformal transformations of the affine 
parameter, the geodesies of the metric g/^i, given by (18) with affine parameter A are the same as 
those of the metric g^i^ = a~'^g^y with affine parameter t] such that d\ = a'^dr). The trajectory 
of the photon coming from the direction given by the Minkowski-nuH vector = {l,n) with 
= 1 wiU be perturbed in such a way that we can write x'^{r]) = n^rj + Sx^, where the second 
term corresponds to the contribution from hij which is of second order. Then, assuming that 
the observer velocity is of first order, the momentum of the photon to second order is: 

(22) 




Inserting this expression in (20) we obtain: 

£^ — \l + -v -v-n-\ ; — . (23) 

a \ 2 dr) J 

For the g^j^i, metric, the second order of the zero component of the geodesic equations in terms 
of the metric perturbation reduces to: 

d^5x^ Idhii • , 

which can be easily integrated to get: 

— — = — hi-jTi'n^. (25) 
dr) 2 ^ ^ 

Then, the energy of the photon results finally: 

E , 1.2 ^ 1 



£=-[l + -v'-v-n--hijn'nn. (26) 

Thus, by using this formula for the energy of the photon in equation (19) and expanding up to 
second order we obtain the following expression for the temperature fluctuations: 

'^^ ^lL-(^<iec-n)(i;.n)|L-^/^iinV|L- (27) 
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The first term in (27) only contributes to the monopole and the second term is a Doppler effect, 
although notice that since the velocities appearing in v-n\^^^^ are referred to the S* = frame, in 
the case of moving fiuids, the dipole is due to the motion of emitter and observer with respect 
to the cosmic center of mass, as commented above. Finally, the last two terms do contribute 
to the quadrupole. The first of them does not depend on the metric perturbation since it 
comes from the second order expansion of the denominator in (19). As we will show below, 
this term is expected to be smaller than the last one. Therefore, the dominant contribution for 
the quadrupole is given by the following expression: 

ST 1 

= -2(^(^0) - hij{adec))n'n^ . (28) 
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This formula shows that we only need to know hij in order to calculate the quadrupole and, 
besides, this term does not depend on the observer velocity. Moreover, it is easy to see from 
(12) that the solution for hij is the following: 



ha = 6 



Ha 'Pa J y^ai ^aj ^^a 




dd dd 



(29) 



where a* is the value of the scale factor at the time at which we specify the initial conditions 
for hij. Notice that the quadrupole does not depend on hij{a*), but only on the initial value 
of the derivatives. For simplicity we will assume that the metric anisotropies are generated by 
the fluids motion and therefore we consider a purely isotropic universe for a < a*, i.e., we take 

hij{a^) = hij{a^) = 0. 

As we can see in the last expression, the quadrupole depends on both the zeroth order 
densities of the fluids and the first order of the velocities with respect to the center of mass. 
These quantities can be obtained from the conservation equations (13) and (15). The zeroth 
order equations have the well-known solutions for the densities p^^^ = poaQ~^^'"""'"^\ Poa being 
the densities today. On the other hand, this expression for the zeroth order densities allows us 
to obtain the solutions for the velocities: 



where voa arc the velocities of the fluids today. We see that, to flrst order, each fluid will move 
along a fixed direction. Well inside the radiation dominated era, baryons and dark matter 
particles were coupled to radiation which, being the dominant component, will drag matter in 
such a way that the three fluids velocities are the same. Prom decoupling on, each fluid velocity 
will start evolving according to (30), i.e. matter will start reducing its velocity with respect 
to the cosmic center of mass frame, whereas radiation keeps moving at a constant velocity. 
Since dark matter is expected to decouple before baryons do, and both velocities scale in the 
same way, the two matter fluids are expected to be moving with constant relative velocity after 
recombination. Finally since the initial velocity of the three fluids (radiation, baryons and dark 
matter) were the same, the direction of their velocities will also agree after recombination. On 
the other hand, the gauge condition 5" = in (6) yields the constraint: 



so we can conclude that dark energy should also move in the same direction as the rest of 
fluids in this frame. That way, each energy-momentum tensor (and therefore the total one) 
will have axisymmetry so the metric will also be axisymmetric. This means that the tensor 
perturbation hij is diagonal. In fact, if we choose the velocities along the 2;-axis, the tensor 
perturbation given by (29) will be proportional to diag{— 1, —1,2). Taking into account the 
previous discussion, the flnal expression for the quadrupole results: 



^ 1 -> 



.3Wa — l 



(30) 




(31) 



a 




(32) 
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where 6 is the angle formed by the observation direction and the velocities of the fluids, and 
h{a) = J2a ha{ci) with 



Ja* Ci J a* 



da 



(0) 



da 



(33) 



The function of 9 appearing in (32) is proportional to the spherical harmonic so we can 
express the quadrupole fluctuation as: 



5Tq 2 /7f 



T 



2 [W 

2^5(^0-^^60)^20- 



(34) 



It is usual to introduce the power spectrum of the temperature fluctuations of CMB as: 



1 £(£ + 1) 
T ~^27r 2£ + l t^''''^ 



(35) 



where a^m are the coefficients of the expansion in spherical harmonics. Moreover, the quadrupole 
is usually deflned as: 



Q 



T \5tt 



m=-2 



which, in our case, reduces to: 



5V3 



(36) 



(37) 



The quadrupole given by (37) is due to the anisotropy of the space-time background (that is 
why the index A is introduced), but we have to add the standard isotropic fluctuation produced 
during inflation. Then, if we assume that the anisotropies are small, the total effect will be the 
linear superposition of both contributions (see [14]): 



STt = STa + 6Tj 



(38) 



and, therefore: 

Notice that, as discussed in [14], there is the possibility that the inflation-produced contribution 
could be strongly biased or anti-biased by the anisotropic background, mainly in the case in 
which anisotropies grew as we go back in time. However as discussed in that reference, this 
is unlikely in general since it would require a correlation between the quantum origin and 
subsequent classical evolution. Moreover, in our case, the background evolution during inflation 
is isotropic and we do not expect any interference effect. 

Following [21], we can easily generalize our results to the case of an arbitrary orientation 
of our frame in which the velocities he along the direction given by {§, 0) . In that case, the 
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coefficients of the expansion are: 



^20 



"21 



*22 



6^5 



[1 + 3 cos 2^] \ho - hdecl , 



-(«2-i)* = -^^ e-^'^sin2^ \ho - hdec\ , 



TT 



= («2-2)* = e-^^* Sin^ e \ho - hdecl ■ 



(40) 



It is easy to show that the anisotropy quadrupole according to (36) is still given by (37) since 
h is scalar under rotations. Now, assuming that the coefficients only differ one from each 
other in a phase factor we can write: 



^20 ~ 



"21 



^22 



3 



(02-2) 



Qi, 



Qi, 



Qi 



(41) 



which is justified because the standard inflation fluctuations are statistically isotropic. Then, 
the total quadrupole can be expressed as: 



(42) 



where / is a function depending on the direction of the velocities {9, 0) and the phase factors 
ttj of the coefficients 03^^, and whose expression is: 



/ 



1 



4v^ 



2^6 2(f) + a^) + 2 cos 9 cos(0 + sin ^ — (1 + 3 cos(2^)) cos cti . (43) 

This function takes values such that: 

l/l</ = 



39 + 6^6 + ^6-1 



4a/5 



(44) 



Since the values of the phases ctj are random, the total quadrupole lies between and Q^, 
being: 



Ql^Ql + Qj± 2fQAQi 



(45) 



with / the maximum of /. 



The observed quadrupole from WMAP [3] is given by {ST)l^, = 236^^1? fxK^ at the 
Vo C.L. or (5T)2^, = 236l?i^ A«K^ at the 95% C.L. These results define the correspond- 
ing 68% C.L or 95% C.L. intervals for the measured temperature fluctuations that we denote: 



((5T)^j„, (5T)^^^). For the theoretical quadrupole temperature interval ((^T)^, {ST)\), ob- 
tained from (45), to be compatible with observations, we therefore require (5T)^^^ > (5T)^ 
and (ST)^-^ < (5T)^. Using (45) these two conditions impose limits on {STaY once the value 
of {5T)j is fixed. 



9 



Let us first assume that inflation alone is able to account for the observed quadrupole, i.e., 
((5T)j ~ then the first condition ((5T)^„^ > (5T)^ is automatically satisfied, because 

the minimum of (ST)'^, as a function of 5Ta, is {ST)^ which is larger than {5T)1^^^. Therefore, 
we obtain bounds on STa just from the second condition above, which are given by: 

O/zK^ < (5Ta)^ < 1861 /xK^ 68%C.L. 

O/zK^ < ((5T^)2 < 5909 /^K^ 95% C.L. (46) 

However, it is well-known that the predictions of standard infiation, calculated from an 
almost fiat spectrum of density perturbations, is larger than the central value of the measured 
quadrupole, in particular: (5T)^ ~ 1252/iK^. In such a case the anisotropic contribution could 
help reducing the value of the quadrupole for certain values of the phases and fiuid velocities. 
Once again the first condition is automatically satisfied, and the second condition yields: 

54 /iK^ < {dTA^ < 3857 /xK^ 68% C.L. 
/xK^ < {STa^ < 9256 /xK^ 95% C.L. (47) 

Notice that the 95% confidence interval includes the standard prediction from infiation and 
for that reason the lower limit vanishes in that case in (47). According to these results, for 
certain orientations of the velocities and the values of the phase factors, Qa could lower the 
value of quadrupole and make it compatible with the observed one even at the la level. We 
will compare these limits with the predictions from several models, but before that we need to 
extend our calculations beyond the perturbative regime. 

4 Fast-moving fluids: exact equations 

In the previous Section we have studied the problem of obtaining the quadrupole produced by 
the fact that dark energy does not share a common rest frame with matter and radiation. To 
that end, we have used cosmological perturbation theory to compute the metric perturbations 
by means of the simple formula (29), valid when the velocities are small. Such formula could 
also be reasonably useful for high initial velocities provided they drop in time and rapidly reach 
the perturbative regime. If we look at Eq. (30) we conclude that this condition is satisfied if 
Wa < |. Besides, if = |, as in radiation case, the velocity is constant so we just need to have 
a small initial velocity. However, some models have been proposed in which the total energy 
density of the universe could contain in certain epochs a non-negligible contribution of fiuids 
with equation of state such that Wa > \- This is for instance the case of stiff-fiuid cosmologies 
or some tracking dark energy models. In those cases the velocities grow in time, perturbation 
theory will eventually breaks down at some point and it becomes necessary to solve the exact 
problem. 

In order to simplify the equations in this case, we shall change the gauge used in the previous 
Section by one in which ^fQo = 1, keeping the condition goi = 0. Moreover, we shall also assume 
that the fiuids are moving along the z-axis with no rotation which means that axisymmetry 
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holds. The most general metric having that symmetry with this gauge choice can be written 
as follows: 



ds"^ = df - alidx'^ + dy'^) - a\dz^ (48) 
On the other hand, the energy-momentum tensor for each fluid reads: 

^aO = 7a(Pa+Pa) -Pa, 

Tai = -7a (pa + Pa) afu^Ji^, 

= -ll{Pa+Pa)a\vlJ^^5,,-pJ]. (49) 

Note that the velocities appearing in these expressions are no longer the same as those of the 
previous Section since, here, we have defined them as derivatives with respect to the time t 
not with respect to r]. However, it is easy to translate these velocities into the others just by 



defining a mean scale factor a = \J a\a\\ because dt ~ adrj and, therefore, ^ ~ which 
is a good approximation in the perturbative regime. Moreover, we have to notice that for an 
appropriate definition of fluid velocity, we have to rescale V = a\\v, so that V"^ < 1. 

Now, it is convenient to introduce the variables Oa defined by cosh 6*0, = 7a so that the 
equations adopt a simpler structure. Velocities are related to 6^ by means of tanh 6'a = ci\\Va- 
With these new variables, the Einstein equations from (48) and (49) take the form: 

Hl + 2H±Hii = 8ttGY, {cosh^ 9a + Wasmh'^ 9a) Pa, (50) 

a 

H± + H\\+Hl + Hf^+H^H\\ = -SnG^pa, (51) 

a 

2H^ + 3Hl = -87rG' ^ (wa cosh^ ^„ + sinh^ (52) 

a 

where ' = ^ and iJ^ = a±/a±, /Jy = a||/a|| are the transverse and longitudinal expansion 
rates respectively. These equations reduce to the Friedmann ones when a± = a\\ and v^a = 0. 
Again, as in the perturbative case, we need some extra equations to close the problem which 
are those given by the independent energy-momentum tensor conservation. These equations 
can be written as follows: 

((Wa - 1) COsh^ 9a - l) H\\ + 2WaH^ 

- ^ ^ —v,a, (53) 



^za — I i\ 1,2/1 

\Wa — 1) cosh 9a — Wa 

(1 + Wa) (H\\ + 2H^ cosh^ 9a) 

Pa = ^TT TJJ. '-Pa. (54) 

[Wa — Ij cosh 9a — Wa 

Besides, one can find the following equations for the evolution of 9a- 

9a = — ^^""^^ sinh 9a cosh 9a- (55) 

[Wa-\) cosh 9a " Wa 
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The spatial geodesic equations for the metric considered are: 

0, 



dX dX 



d'^y 2^ dy ^ 
dX^ ^ dXdX ' 



d'^z dt dz 



0, 



(56) 



where A is an afhne parameter. The first integral for these equations is given by: 



df I Uy Uy, 



dX 



2 ' 2 ' 2 

ai ai ai 



(57) 



being f— (a;, y, z) and the integration constants can be chosen for simplicity in such a way that 
— 1. Moreover, from the condition of null geodesic we get: 



dt _ 

dX ^ a 



(58) 



with n]_ — n^ + and n| = n^. Then, for an observer with velocity — 7(1, v) the energy of 
the photon is: 



and the corresponding temperature fluctuation reads: 



-\ — ^ — n • w 
ai af, 



(59) 



5T 



70^0 



M I II -> -> 



- 1 



T<iec ^dec / n'^ 



(60) 



-'Xdec 



+ TT^ fi ■ Vdec 



where again the indices and dec denote the present and decoupling times respectively. 



5 Model examples 

5.1 Constant equation of state 

The simplest dark energy model we will consider is that corresponding to a fluid with constant 
equation of state wde — —1- Note that in the case wde — —1, i-e. pure cosmological constant, 
dark energy does not contribute to the center of mass velocity (6) and, therefore, the center 
of mass frame agrees with the radiation frame. This means that all the fluids would share a 
common rest frame and no effects on the CMB would be possible. When wde is close to — 1, 
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the velocity of dark energy scales as 



and its energy density is nearly constant (see Fig. 



1). Since dark energy velocity decreases very fast, its contribution to the quadrupole is very 
small. Besides, the velocity of radiation (and therefore that of matter) is determined by the 
initial dark energy velocity and the gauge condition (31) as: 



vr 



1 + WpE ^D E 

1 + wr Q} 



4„-*K 



DE- 



(61) 



where v}-,^ is the initial dark energy velocity. Taking VLj:,e = 0.73, = 8.18 x 10~^, wde = 
—0.97 and a* ~ 10~^, we get -iTr ~ 2 x 10~'^^v}-,^. The value of a* taken corresponds to a 
favorable case since lower values would lead to much lower velocities of radiation and matter. 
Then, even for initial velocities of dark energy close to 1, the velocities of matter and radiation 
are extremely small, which means that the three fluids are very nearly at rest in the cosmic 
center of mass frame. That way, the quadrupole generated in this model is totally negligible. 
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Figure 1: Evolution of densities and velocities in a model with constant equation of state as described 
in the text. Continuous line (blue) for dark energy, dashed-dotted (red) for radiation, dotted (cyan) 
for dark matter and dashed (green) for baryonic matter. On the left baryonic and dark matter are 
added together and plotted in dashed (green). Notice that in this plot dark matter is assumed to 
decouple at z ~ 10^. 



5.2 Scaling models 

Scaling models [22] are those with equation of state such that dark energy mimics the dominant 
component of the universe throughout most of the universe evolution. Thus, dark energy evolves 
as radiation before matter-radiation equality and as matter after that. However, in order to 
explain the accelerated expansion of the universe, dark energy has to exit from that regime and 
join into one with wde < —1/3 at some point. Then, the evolution of the dark energy density 
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is given by: 



PDE 



Pdeo Ot 



3WDE 



Qeqa 



a < 



eq 



Pdeo c^t ^ 



a^q < a < qt 



(62) 



Pdeo o 



a > qt 



where as commented before, ar is the scale factor when dark energy leaves the scaling regime 
and Pdeo is the present value of the dark energy density. 

In the evolution of dark energy velocity, we have to take into account the momentum 
conservation equation given (to first order) by (15). This equation implies that the dark energy 

velocity must be discontinuous at the transition points since the equation of state jumps at 
those times whereas the quantity a'^(l + wde)pdeVde is constant, being pde continuous. With 
this in mind, we get the following evolution for dark energy velocity: 



vde 



V* 



DE 



2 ^eq^ 



WE 



a < a 



"eq 



eq 



< a < Gt 



(63) 



-a 



The discontinuities in the velocity arise because we are considering abrupt changes in the 
equation of state. If these changes were smooth, the results would be essentially unaffected 
since the final values of the velocities would remain those in (63). In Figure 2, we show the 
evolution of the energy densities and velocities for a typical scaling model. 

Wc can see from the previous expression that, in the second transition, the closer tv^E is 
to —1, the more the velocity grows after the transition. The case woe = —1 is not divergent 
because, if that was the case, the conservation equation would become trivial and the velocity 
evolution got from (30) would not make sense anymore. 

In these scaling models, the first transition can be set at the matter-radiation equality 
and the second one must be chosen such that we get the observed dark energy density today. 
Moreover, the initial velocity of radiation, and therefore that of matter, is fixed by the initial 
velocity of dark energy via the gauge condition (31). Because matter is subdominant with 
respect to radiation before equality, the matter contribution in (31) can be neglected and we 
obtain that: 



-ev 



DE 



(64) 



where e = Pde^cl*)/ Pr^cl*) is the initial dark energy density fraction (neglecting the matter 
contribution). Notice that this fraction does not depend on a*, because dark energy scales as 
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Figure 2: Densities and velocities evolution in a scaling model with vj^^ = 0.1 and e = 0.1. As in 
the previous figure, the continuous line (blue) is for dark energy, dashed-dotted (red) for radiation, 
dotted (cyan) for dark matter and dashed (green) for baryonic matter. On the left baryonic and dark 
matter are added together and plotted in dashed (green). Notice that in this plot dark matter is also 
assumed to decouple at z ~ 10^ and = 10~^°. 



radiation in the radiation dominated era. Then, we can obtain a relation between and e just 
by computing that quotient from the known expressions for the energy densities evolutions of 
each fluid. When doing that it results: 



eq 



(65) 



Since we need wj:,e{zt) < —1/3 in order to have accelerated expansion, we see from the 
previous formula that grows as e grows, more precisely if we take wj:,e{z < zt) = —0.97, 
then qt oc e°'^^. Since primordial nucleosynthesis imposes an upper limit on e, we can establish 
also an upper limit on just by setting the maximum value of e on (65). This maximum value 
is emax — 0.2, (see for instance [23]) so we get the constraint ar < 0.41 for this kind of scaling 
models where we have taken VLde = 0.73, VLr = 8.18 x 10~^ and a^q = 

By inserting the corresponding values for the densities and velocities of the four fluids, as 
well as the equation of state considered into (33), we can compute the quadrupole produced by 
the relative motion of the fluids. For our calculations we shall take fls = 0.046, ^Idm = 0.23, 
Qdec = lYoo WDsiz < Zt) = —0.97 and the values given above. Fig. 3 shows the evolution 
of the contribution of each fluid to the metric perturbation ha- We can see that the typical 
behavior is a rapid growth during the radiation era to reach finally a slightly growing regime in 
the matter era (notice that the dependence on a^, is only logarithmic). In spite of the fact that 
the perturbation is 0(f|)|;), the quadrupole is expected to be smaller because h barely grows 
in the epoch since decoupling to today and, as we mentioned in Section 1, the quadrupole is 
essentially given by the growth of the perturbation during that epoch. 
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Figure 3: Evolution of metric perturbations due to each fluid in a scaling model with e = 0.1 and 
f* = 0.1. We can see that the largest contribution comes from dark energy (continuous blue line) and 
matter (dashed green) essentially does not contribute. Radiation is shown with a dashed-dotted (red) 
line. 



The quadrupole produced by scaling models is fixed by two parameters: the initial velocity 
v^^ and the initial energy ratio e of dark energy. It is easy to see from (33) that and, 
therefore the quadrupole, is proportional to f|)|;. Obviously, this dependence is valid just for 
small velocities since when we consider velocities close to the speed of light we have to take into 
account relativistic effects. The dependence of the quadrupole on e can be found to be linear 
for e < 0.07 with a slope 0.44 so we can conclude that the quadrupole is very well approximated 
by the simple expression: 

QA^O.Uev*j,%. (66) 

As commented before, this expression is valid only for small velocities. According to the bounds 
on the quadrupole obtained in (46) and (47), there are allowed regions in the parameter space 
(e, v^^), which from (66) are limited by the curves e = k±/v1 where the constants k± correspond 
to the upper and lower limits on Qa- In Fig. 4 we show these regions obtained numerically 
with the exact equations (see Section 4). As we said above, the second order calculation is a 
good approximation for velocities lower than 0.1. However, when the velocities are large (close 
to 1) values of e < C(10~^) are necessary in order to explain the observed quadrupole. Notice 
once again that these regions have been obtained in the case in which the measured quadrupole 
has two contributions, one coming from inflation and a second contribution coming from the 
fluids motion. 

To end this section we shall show why we can neglect the third term in (27) with respect 
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Figure 4: Exclusion plot in the parameter space (e, v"^^) for a scaling dark energy model. The allowed 
region corresponds to the limits given in (46). The dark (red) strip corresponds to the regions for 
which Qa could explain the observed quadrupole at the 68% C.L. according to (47) 

to that containing the metric perturbation hij. Let us recall that term: 

iv-n)\deci^dec-n). (67) 

The first factor in this expression is nothing but the dipole which is ~ 10"'^. The second factor 
contains the velocity of the observer at decoupling time which coincides with matter velocity 
(and therefore with that of radiation) at that moment. Then, if we recall the relation (64) 
between radiation and dark energy velocities, we find that this term is ~ 10~^ev}^^. On the 
other hand, the last term in (27) is ~ ef|)|; as we have just seen above. Hence, if we call 
and Qh to the last two terms in (27) respectively we have that Qh ~ lO^v'^^Qv and we see 
that will be larger than Qh only for v'^^ < 10~^. However, in such a case the contribution 
to the quadrupole is < 10~^e which is negligible. 

5.3 Tracking models 

In this Section we would like to comment on the difficulties which can appear in certain dark 
energy models when we consider perturbations in the fluids velocities. In general, any model 
with a stiff stage in which its equation of state satisfles w > |, would be unstable with respect to 
velocity perturbations according to (30). This could be the case of certain tracking models [24]. 
These are models in which the energy density of dark energy follows a common evolutionary 
track for a wide range of initial conditions. This attractor behavior makes this kind of models an 
interesting alternative to a cosmological constant since they alleviate the so called coincidence 
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problem. Unlike scaling models, in this case dark energy does not necessarily mimics the 
dominant component. In the model proposed in [24] the equation of state is initially close to 
1, then it changes to —1 and, finally, it oscillates around —0.2. Fig. 5 shows a typical behavior 




Figure 5: Densities and velocities evolution in a typical tracking model with an initial equation of 
state wde = 0.9, which changes to wde = —1 and then to wde ~ —0.2. The continuous line (blue) 
is for dark energy in the model with moving dark energy, whereas the blue dotted line is for static 
dark energy, dashed-dotted (red) for radiation, dashed (green) for baryonic and dark matter. We see 
that when Vde reaches 1, the corresponding density begins falling too fast to be able to recover the 
present value for JIds- 

when w > first the velocity perturbation (defined as y = a\\v) grows according to (30) and 
asymptotically approaches 1. We can understand this from the exact conservation equations 
(53) by taking the ultrarelativistic limit 6^1. This yields the solutions: 

V = voaj^, 

p = po«± . (68) 

This means that \^ ^ 1 is a solution of the equations. In addition, the energy density for w > 1/3 
falls very fast with the expansion when compared with the usual behaviour p = poa~^^^~^'^\ In 
addition, ■y'^p, which is the quantity that contributes to the Hubble rate in (50), decays as 
{a\\a±)~'^, once the fiuid reaches the ultra-relativistic regime, regardless the value of w. 

In the limiting case of stiff fluids with w = 1, it is possible to obtain exact solutions. Thus 
the velocity perturbation and the energy density are: 

V = Voal, 

g-4 J cosh^ edt 

P = Po 2 ■ (69) 

a| 

Thus, the velocity of the fluid grows as a\ until it reaches the speed of light in a flnite time 
and the density falls to zero at the same time because 9 becomes inflnity at that moment. From 
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that time on, the fluid will keep moving at the speed of light with vanishing energy density. 
Notice that 7 becomes infinity in such a way that the momentum is conserved i.e. the mathcing 
between the two regimes must be taken so that is finite and continuous. In any case, we 
see that when the velocity is high enough the density falls to zero and, the closer is w to one, 
the faster is the fall of the density, so that we cannot recover the present value for the dark 
energy density. Notice that this general behaviour is independent of the value of the initial 
velocity and, accordingly, even in models in which all the fluids are initially at rest, a small 
perturbation in the velocity could change dramatically the flnal values of the densities, unless 
flne tunings of the transition redshifts are introduced. 



5.4 Null dark energy 

In this Section we shall study the case in which dark energy behaves as a null fluid, whose 
energy-momentum tensor reads: 

^{Pn+PnTI" -PNg^" (70) 

with a null vector given, in the Bianchi type I metric (48), by P = (1, 0, 0, a||"^). For this 
kind of fluid, the conservation of energy and momentum can be expressed as follows: 

= p^, (71) 
= {pn+Pn) + 2{H\\+H^){pn+Vn)- (72) 

These equations imply that the pressure is constant and that (pat + 'Pn) scales as (a||aj_)~^, so 
that the energy density is given by pat = pm{ci\\(^i.) '^ ~Pm where patq and pm are constants of 
integration. Since the anisotropy is expected to be small, the energy density of this fluid behaves 
as radiation during the early epoch and as a cosmological constant with energy density —pNo at 
late times. Now, if we require pat to be positive at all times, we conclude that the pressure must 
be negative, as corresponds to a cosmological constant. Notice that this is a general result for 
any null fluid whose energy-momentum tensor is given by (70). The transition between both 
regimes can be easily calculated and it is given by ay — {~^^Y^^- Since pm — — 0.73pco, 
where Pco is the critical density today, we have that a-r = O.le^/^ where e = ^ is the ratio of 
dark energy density with respect to radiation which is almost constant. Besides, this ratio is 
also the initial contribution of dark energy to the total energy density which has an upper limit 
imposed by primordial nucleosynthesis. Taking once again Cmax ^ 0.2, we have an upper limit 
on the transition given by: ot < 2 x 10~^. 

The exact Einstein equations in this case are: 

Hl + 2H^H\\ = 87rG'^(cosh2^„ + 'u;„sinh^^„)p„ + 87rG'pjv, (73) 

a 

H^ + H^^ + Hl + Hf^+H^H^^ = -87rGJ2Pa-8TTGpN: (74) 

a; 

2H± + 3Hl = -SnG ^ (wc, cosh^ 9^, + sinh^ Oa) Pa - 8t^G{pn + ^p^) , 

a 

(75) 
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where now a = B, DM, R. Moreover, we still have the gauge condition S = which yields the 
following constraint: 

J2 + Pa)Va + {pN + Pn) = 0. (76) 

a 

In the radiation-dominated era we can neglect the contribution of matter (dark matter and 
baryons) to the latter sum, so we get: 

2 Pn +Pn f-,-,-. 

Pr+Pr 

Since, in that epoch, dark energy must be subdominant with respect to radiation, the quotient 
on the RHS is small and, therefore, the velocity of radiation is also small. This allows us to 
consider the perturbative regimen in the velocities (except, obviously for the null fluid). 

Therefore, if we assume that the anisotropy generated is small we can set the following form 
for a\\ and aj^: 

a_L = a(l + 5_L), 

a\\ = a(l + 5||). (78) 

With this ansatz it is easy to see that h = 2{d\\—S±). Then, inserting (78) in (75) and expanding 
up to first order in 5's and Va we can get the following equation for h: 

|(a»f)^2<.3(,„ + p„). (79) 

This equation can be easily solved by means of two direct integrations and its solution can be 
expressed as follows: 



a» Lb 



/ a^ipN+PN) 

J a* 



dd 



VT^aPa 



dd , , 

(80) 



V^a Pa 



In principle, the problem is not solved yet since pn + Pn depends on and 6±. However, we 
can consider the lowest order in this quantity, i.e., pN + Pn = PnocT^ to obtain the dominant 
contribution to the quadrupole. This is justified because ^ is of the same order as vr as we 
can see from (77). We have to note that, to this order, the quadrupole depends just on the 
null fluid because the first contribution to the anisotropy due to the rest of fiuids is of second 
order in the velocities, whereas the null fiuid contributes to first order. In Fig. 6 we plot the 
evolution of the fiuids densities and h function for a null fiuid with e = 5 x 10"^. 

In this model we only have one free parameter: p^o or, equivalently, e = so we can get 
bounds on e just from (46) and (47). Besides, the quadrupole is linear in e, as we see looking 
at (80), more precisely we have that the quadrupole is given by: Qa — 2.58 e. Now again the 
contribution from the term Q„ given by (67) is negligible compared with Qa since ~ 10~^e. 

This expression is nearly independent of a* because the quadrupole depends on the difference 
ho—hdec which is not very sensitive to the time at which we set the initial conditions. Comparing 
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Figure 6: Left: densities evolution for a null fluid with e = 5 x 10^^. Matter is plotted with dashed 
(green) line, dotted-dashed (red) for radiation and continuous (blue) line for dark energy. We see 
that the null fluid behaves as in a scaling model except during the matter dominated epoch. Right: 
evolution of h showing that the anisotropy grows up to a maximum value where it keeps almost 
constant. 

the expression obtained for the quadrupole with the previous bounds, we get that the allowed 
region in (46) corresponds to: 

e < 6.1 X 10"^ 68% C.L. 

e < 1.1 X 10"^ 95% C.L. (81) 

whereas for: 

1 X 10"^ < e < 8.8 X 10"^ 68% C.L. 

< e < 1.4 X 10"^ 95% C.L. (82) 

the null fluid could make the predicted quadrupole to agree with observations, as shown in 
(47). 

6 Conclusions and discussion 

In this paper we have studied homogeneous models of dark energy in which the rest frame of 
the different fluids can differ from each other. We have considered the evolution of slow-moving 
and fast-moving fluids and shown that, starting from an initially isotropic universe, the fluids 
motions can generate an anisotropic expansion in which the anisotropy degree typically grows 
in time. Such anisotropics are shown to contribute to the CMB dipole and quadrupole only, 
but not to higher multipoles. We apply those results to some dark energy models and flnd that 
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in models with constant equation of state, even for initial velocities of dark energy close to the 
speed of light, throughout the matter era all the fluids would practically share a common rest 
frame and therefore no effect on the quadrupole is expected. However, in the case of scaling 
models it is shown that the anisotropy grows during the radiation era and reaches an almost 
constant value during matter domination. The effect on the CMB quadrupole can be relevant 
and bounds on the velocity and initial fraction of dark energy can be found. We also find that 
for models with an initial stage in which the equation of state is stiffer than radiation, as for 
instance in some tracking models, the velocity approaches the speed of light whereas the energy 
density decays faster than in the case in which dark energy is at rest with respect to matter 
and radiation. This fact spoils the predictions of those models for the density parameters at 
late times. Finally we have considered also fluids moving at the speed of light and found that 
generically they behave as a cosmological constant at late time, provided their energy density 
is positive at all times, whereas they act as radiation at early times. The contribution to the 
quadrupole is also used to set limits on the relative contribution of dark energy in the radiation 
dominated era. 

The models presented in this paper with moving fluids lead to Bianchi I anisotropic metrics 

of the type recently studied in [21]. However unlike that model with decaying anisotropics 
generated by the presence of magnetic fields in the early universe, in our work the motion of 
dark energy supports the anisotropics which could have a non-negligible value today. 

Finally concerning the potential effects on the CMB polarization, it is well known [25] that 
Bianchi models give rise to polarized radiation through Thomson scattering at decoupling time. 
This is also the case during the reionization period and therefore a potential contribution from 
dark energy motion is expected also in the polarization signal at large scales. Work is in progress 
in this direction. 
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